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Nonlinear Asymptotic Stability of the Lane-Emden Solutions for 
the Viscous Gaseous Star Problem with Degenerate Density 

Dependent Viscosities 

Tao Luo, Zhouping Xin, Huihui Zeng 
Abstract 

The nonlinear asymptotic stability of Lane-Emden solutions is proved in this paper 
for spherically symmetric motions of viscous gaseous stars with the density dependent 
shear and bulk viscosities which vanish at the vacuum, when the adiabatic exponent 
7 lies in the stability regime (4/3,2), by establishing the global-in-time regularity 
uniformly up to the vacuum boundary for the vacuum free boundary problem of the 
compressible Navier-Stokes-Poisson systems with spherical symmetry, which ensures 
the global existence of strong solutions capturing the precise physical behavior that 
the sound speed is C^/^-Holder continuous across the vacuum boundary, the large time 
asymptotic uniform convergence of the evolving vacuum boundary, density and velocity 
to those of Lane-Emden solutions with detailed convergence rates, and the detailed 
large time behavior of solutions near the vacuum boundary. The results obtained in 
this paper extend those in [30] of the authors for the constant viscosities to the case of 
density dependent viscosities which are degenerate at vacuum states. 


1 Introduction 


The evolution of a viscous gaseous star in three 
the interface between the gas and vacuum can 
problem 

' pt + div(pu) = 0 

(pu)t -k div(pu (g) u) -k div6 = -pVx^' 

p > 0 

p = 0 and ©n = 0 
V{T{t)) = u ■ n, 

. {P, u) = (po,uo) 

Here (x, f) G x [0, cxo), p, u, & and T denote, respectively, the space and time variable, 
density, velocity, stress tensor and gravitational potential; Q{t) C T{t), V(r(f)) and n 
represent, respectively, the changing volume occupied by a fluid at time t, moving interface 
of fluids and vacuum states, normal velocity of r(f) and exterior unit normal vector to T{t). 
The gravitational potential is given by 

T(x, f) = —G [ dy, satisfying AT = dvrGp in Q{t) 

Jn{t) 1^ “ y| 

with the gravitational constant G taken to be unity for convenience. The stress tensor takes 
the form: 


spatial dimensions with the boundary being 
be modeled by the following free boundary 


in 

n{t), 


in 

Q{t), 


in 

n{t), 


on 

Tit): 

= d^{t), 

on 

Q : = 

H(0). 


e 


p/3 — Ai ( Vu -k Vu* 


-(divu)J3 


A2(divu)J3, 


1 





where I 3 is the 3x3 identical matrix, p is the pressure of the gas, Ai is the shear viscosity, 
A 2 is the bulk viscosity, and Vu* denotes the transpose of Vu. We consider the polytropic 
gases for which the equation of state is given by 


P = P{P) = 

where > 0 is a constant set to be unity for convenience, 7 > 1 is the adiabatic exponent. 

For a non-rotating gaseous star, the stable equilibrium configurations, which minimize the 
energy among all possible configurations (cf. [21]), are spherically symmetric, called Lane- 
Emden solutions. Therefore, spherically symmetric motions for dynamical problems are 
important to consider. When the viscosities are positive constants, the global-in-time spher¬ 
ically symmetric solution to the free boundary problem fll.ip and its nonlinear asymptotic 
stability toward the Lane-Emden solution were proved in [30] for 4/3 < 7 < 2 (the stable 
index), by establishing the global-in-time regularity uniformly up to the vacuum bound¬ 
ary of solutions capturing an interesting behavior called the physical vacuum which states 
that the sound speed c = \/p'{p) is C^/^-Holder continuous near the vacuum boundary (cf. 
[21131 Ha [I6112311251 Sg), as long as the initial datum is a suitably small perturbation of 
the Lane-Emden solution with the same total mass. The large time asymptotic convergence 
of the global strong solution, in particular, the convergence of the vacuum boundary and 
the uniform convergence of the density, to those of the Lane-Emden solution with detailed 
convergence rates as the time goes to infinity are given in [30] when the viscosities are con¬ 
stant. The aim of this work is to extend those results to the more physically reasonable case 
that the viscosities depend on the density and vanish at vacuum states. The degeneracy of 
the viscosities near vacuum is one of the main difficulties which makes such a generalization 
highly nontrivial. We assume that the viscosities Ai and A 2 are solely functions of the density 
in this paper. For simplicity, we set 

Ai = i^ip^ and A 2 = v' 2 p^, where z/i, z /2 and 9 are positive constants. ( 1 - 2 ) 

One may check easily that the theorems proved in this article apply to the case that Ai and 
A 2 are positive for p > 0 and have the behavior of p^ as p —)■ 0 . 

In the spherically symmetric setting, that is, Q{t) is a ball with the changing radius R{t), 

p(pi.,t) = p{r,t) and u(x, t) = u{rR)^/r with r = |x| G {0,R(t)) ; 
system fll.ll) can then be reduced to 




ir‘^P)t + ir^pu)r = 0 


in (0, R{t )), 



in (0, R{t )), 

in [0, R{t )), 
for r = R{t), 


(1.3) 


p > 0 
p = 0 

R{t) = u{R{t),t) with i?(0) = Rq, u{0,t) = 0, 
. {p,u) = (po,Mo) 


on (0, Rq). 


2 












The initial domain is taken to be a ball {0 < r < -Rq}- And the initial density is supposed 
to satisfy the following condition: 

Po{r) >0 for 0 < r < Po(-Ro) = 0 and — oo < ^ ^ ^ ~ 


So, 

^ Rq — r as r close to Ro, (1.5) 

that is, the initial sound speed is continuous across the vacuum boundary, which 

is called the physical vacuum for the compressible inviscid flows (cf. [H |31 [T31 [251 112 ] )• 

The behavior fll.dp near the vacuum boundary captures a very interesting feature of 
the stationary solution of fll.Sp . {p,u) = (p, 0), the Lane-Emden solution (cf. [H |22]). 
corresponding to a non-rotating gaseous sphere in hydrostatic equilibrium. Here p solves 

dr{p^) + ^Tir~‘^p f p{s)s^ds = 0. ( 1 . 6 ) 

Jo 

The solutions to fll. 6 p can be characterized by the values of 7 (cf. [22]) for given finite total 
mass M > 0. If 7 G (6/5,2), there exists at least one compactly supported solution, and 
every solution is compactly supported and unique for 7 e (4/3,2). If 7 = 6/5, the unique 
solution admits an explicit expression, and it has inhnite support. On the other hand, for 
7 G (1,6/5), there are no solutions with finite total mass. For 7 > 6/5, let R be the radius 
of the stationary star giving by the Lane-Emden solution, then it holds (cf. [22l [31]) 

p'^“^(r) R — r as r close to R. (1.7) 

In both astrophysics and the theory of nonlinear PDFs, the problem of nonlinear asymp¬ 
totic stability of Lane-Emden solutions is of fundamental importance. As mentioned in [30] . 
the key to this is to establish the the global-in-time regularity of higher-order derivatives of 
solutions uniformly up to the vacuum boundary, which has been challenging due to the high 
degeneracy of system (ll.3p caused by the behavior (II.5p near the vacuum boundary, even for 
the local-in-time existence theory in both the inviscid and viscous compressible flows. Indeed, 
the local-in-time well-posedness of smooth solutions to free boundary problems with physical 
vacuum was only established recently for compressible inviscid flows (cf. [ 2 ] [S] [TS] [16] , and 
[ 22 ] which proved a local-in-time well-posedness theory in a new functional space for the 
three-dimensional compressible Euler-Poisson equations in spherically symmetric motions). 
For the vacuum free boundary problem fll.3p of the compressible viscous flows featuring the 
behavior fll.5p near the vacuum boundary, a local-in-time well-posedness theory of strong 
solutions was established in [H] and [7], for the constant and density dependent viscosities, 
respectively. In [HU], the authors succeeded in establishing such a global-in-time higher-order 
regularity for the constant viscosities. Compared with the case of constant viscosities, the 
degeneracy of viscosities at vacuum states makes this a much more challenging task. In 
addition, the term 4:{Xi)j.u/r = Ai/i{p^)ru/r in the momentum equation, which does not ap¬ 
pear in the constant viscosity case, causes serious difficulties in the analysis for the solutions 
with the behavior fll.5p . Indeed, near the vacuum boundary, r = R{t), this term involves 
9{R{t) — which is unbounded for 0 < 6 ^ < 7 — 1, where 9 is the constant in fll. 2 p . 
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This difficulty appears even in the study of local existence of strong solutions. In fact, the 
local existence result of strong solutions to fll.3p in [7] only holds for 6 = 1, which avoids 
the unboundedness of Aui{p^)ru/r for 0 < 6 ^ < 7 — 1 when 7 < 2. In the present work, 
special care is taken to deal with this term near the vacuum boundary by constructing suit¬ 
able weights to capture the behavior of velocity near the vacuum boundary to resolve this 
difficulty. In addition to this, we also rehne the arguments in [30] substantially in obtaining 
the L°°-bound for the first derivative of velocity, for which the weighted L^-estimates of the 
second derivative are used in [30], by establishing a pointwise estimate for the first derivative 
of velocity near the vacuum boundary only involving the weighted L^-estimates of the first 
derivative. 

When 3/4 < 7 < 2, we extend in this article the nonlinear asymptotic stability results 
in ED] for the constant viscosities to the case of the density dependent viscosities ( 1 . 2 ) with 
0 < 6 * < 7 / 2 , by proving the existence of a unique global-in-time strong solution to fll.31) 
and establishing the global-in-time regularity uniformly up to the vacuum boundary, which 
ensures the large time asymptotic uniform convergence of the evolving vacuum boundary, 
density and velocity to those of the Lane-Emden solution with detailed convergence rates, 
and detailed large time behaviors of solutions near the vacuum boundary. In particular, we 
show that every spherical surface moving with the fluid converges to the sphere enclosing 
the same mass inside the domain of the Lane-Emden solution with a uniform convergence 
rate. Therefore, the convergence of the vacuum boundary r = R{t) as f —?■ 00 to that of 
the Lane-Emden solution r = R holds as a consequence. This also shows that the large 
time asymptotic states for the vacuum free boundary problem fll.311 are determined by the 
initial mass distribution and the total mass. The results obtained in the present work are 
among few results of global-in-time strong solutions to vacuum free boundary problems of 
compressible fluids capturing the singular behavior of fll.51) . 

Extensive works have been done on the studies of the Euler-Poisson and Navier-Stokes- 
Poisson equations with vacuum, especially in recent years. One may find the study of 
the stability problem of gaseous stars in astrophysics literatures (cf. [H HU [IS])- The 
linear stability of Lane-Emden solutions was studied in [22|- By assuming the existence of 
global solutions of the Cauchy problem for the three-dimensional compressible Euler-Poisson 
equations which has been a major challenge in the theory of fluid dynamics equations, a 
conditional nonlinear Lyapunov type stability theory of stationary solutions for 7 > 4/3 
was established in [SB] using a variational approach (the same type of nonlinear stability 
results for rotating stars were given by [2H1127I). Those nonlinear stability results are in the 
framework of initial value problems in the entire and involve only Lyapunov functionals 
which are essentially equivalent to L^-norms of the difference of solutions, where the vacuum 
boundary cannot be traced. In the framework of free boundary problems for the Euler- 
Poisson and Navier-Stokes-Poisson equations, the nonlinear dynamical instability of Lane- 
Emden solutions for 7 G (6/5,4/3), was proved in [TB] and [T7|, respectively. In inviscid 
flows, a nonlinear instability for 7 = 6/5 was proved in [12]; an instability was identified for 
7 = 4/3 in [1] that a small perturbation can cause part of the mass to go off to inhnity. 

It should be noted that the existence of global weak solutions was proved in [10] for the 
initial boundary value problem reduced from the vacuum free boundary problem fll.3p after 
using the Lagrangian mass coordinates, under the constraint that 7 > 4/3, 9 E ( 0,7 — 1) fl 
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( 0 , 7 / 2 ], and 


2(8 + 4a-a2) 2 8^5 + 4a - V 2 2(8 + 4a - a^) 2 8^5 + 4a - a^ 

4 — 4a + a^ 3 4 — 4a + a^ ^ z/i ^ 4 — 4a + a^ 3 4 — 4a + a^ 

with a G (—1,1) being a constant. In contrast to the strong stability result shown here, for 
the global weak solutions obtained in [ 10 ], only the uniform convergence of the velocity is 
proved, due to the lack of regularity near the vacuum boundary, and the uniform convergence 
of the density, in particular, the convergence of the vacuum boundary which is the most 
interesting part in the study of asymptotic behavior of the free boundary problem, are 
missing. Furthermore, our nonlinear asymptotic stability result holds for 4/3 < 7 < 2, 
0 < 0 < 7 / 2 , ui > 0 and 1^2 > 0 , without the restrictions on 0 < 7 — 1 and the ratio of ^’ 2 / 1^1 
as in [To]. It should be noted that 7 — 1 < 7/2 for 7 < 2 . 

We conclude the introduction by reviewing some previous works on viscous flows. It 
should be noted that there are also other prior results on free boundary problems involving 
vacuum for the compressible Navier-Stokes equations besides the ones aforementioned. One 
may refer to [33l [SS] [28l [9l [181 El SSI IHl UHl |15] and references therein for the one-dimensional 
motions concerning global weak solutions. For the spherically symmetric motions, global 
existence and stability of weak solutions were obtained in [311 ES] for gases surrounding a 
solid ball (a hard core), restricted to cut-off domains excluding a neighborhood of the origin; 
a global existence of weak solutions containing the origin was established in HH for which 
the density does not vanish on the boundary and the viscosities are density dependent. 
The readers may refer to mm for the local-in-time well-posedness results and [iQl for 
linearized stability results of stationary solutions for a class of free boundary problems of 
the compressible Navier-Stokes-Poisson equations away from vacuum states. 

2 Main Results 

First, we recall some properties for Lane-Emden solutions. For 7 G (4/3,2), it is known 
that for any given hnite positive total mass, there exists a unique solution to equation fll.bp 
whose support is compact (cf. [22]). Without abusing notations and for convenience, we use 
X as the variable in the study of Lane-Emden solutions. That means, for any M G (0, cxo), 
there exists a unique function p{x) such that 

po := p(0) > 0, p{x) >0 for a; G (O, R) , P {R) =0, M= ATTp{s)s‘^ds-, (2.1) 

Jo 

—00 < px < 0 for a: G (0, R) and p{x) < po for x G (O, R) ; (2.2) 

(fR)^ =—x(j)p, where (j) := x~^ f 47 rp(s)s^(is G [M/R^, 47 rpo/ 3 ] ; (2.3) 

ao 

for a certain hnite positive constant R (indeed, R is determined by M and 7 ). Note that 

(rh = ( 2 . 4 ) 
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It then follows from fl2.ip and fl2.3p that p satishes the physical vacuum condition, i.e., 


p'^ ^ (x) ~ i? — X as X close to R. 


Indeed, there exists a constant C depending on M and 7 such that 
C-^ {R-x)< p^-\x) <C{R-x), X e (0, i?) . 


(2.5) 


We adopt a particle trajectory Lagrangian formulation for fll.Sp as follows. Let x be the 
reference variable and define the Lagrangian variable r{x,t) by 

rt{x,t) = u{r{x,t),t) for f >0 and r(x, 0 ) = ro(x), xe/;=( 0 ,.R). 

Here ro(x) is the initial position which maps I —?■ [0, i?o] satisfying 

rroix) rx 


/ po{s)s^ds = / p{s)s^ds, X G /, 

Jo Jo 

SO that 

Po(''o(a:))ro(^)''o(^) = x € 1. 


( 2 , 6 ) 


(2,7) 


(Indeed, fl2.6p means that the initial mass in the ball with the radius ro(x) is the same as 
that of the Lane-Emden solution in the ball with the radius x. Then smoothness of ro(x) at 
X = .R is equivalent to that the initial density po has the same behavior near Rq as that of 
p near R.) The choice of tq can be described by 


ro(x) = -0 ^(^(x)), 0 < X < i?; 
where ^ and ijj are one-to-one mappings, dehned by 


( 2 . 8 ) 


^ : {0,R) —)■ (0, M) : X 1 -^ / s‘^p{s)ds and '0 • (0, ^ (0,M) : 2 ; i—)■ / s‘^po{s)ds. 

Jo Jo 

Moreover ro(x) is an increasing function and the initial total mass has to be the same as 
that for p, that is. 


cRo 


4:7rpo{s)s‘^ds = 


i'ro{R) 


r-i? 


4:7rpo{s)s^ds = / 4:7ip{s)s^ds = M, 


(2.9) 


to ensure that tq is a diffeomorphism from I to [0, i?o]- It follows from fll.Sp ^. that 

pr{x^t) prQ{x) 

/ p{s,t)s^ds = / po{s)s'^ds, X G/. 

Jo Jo 

Set 

f{x,t)=p{r{x,t),t) and v{x,t) = u{r{x,t),t). 


( 2 . 10 ) 
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Then the Lagrangian version of system fll.3l) ^ ^ can be written on the reference domain J as 

(rV)i + rV- = 0, 


/■ro(a:) i 

fvt + ^ + 47r/r“^ / po{s)s‘^ds = — 


-Ai + A2 




( 2 . 11 ) 


-4(Ai), 


TxT 


Solving fl2.1ip ^ gives that 

f{x,t)r^{x,t)r^{x,t) = poiro{x))rl{x)ro^{x), x e I. 

Therefore, 




x^p{x) 


for X E I, 


r'^{x,t)r^{x,t) 

due to fl2.7p . By using fl2.3p . the free boundary problem fll.3p is then reduced to the following 
initial boundary value problem on a fixed interval J; 


'X\2 

P[Z) + 


2 — \ 7 ' 
xp^ 

r 2 J. 


X 


+ ^P ^7Tpy dy = V{x,t) in /x(0,T], 


n(0, t) = 0 

^ (r, n)(x,0) = (ro(a;), Mo7o7))) 


on (0,T], 
on I X {t = 0}, 


( 2 . 12 ) 


where 


V = z/ 


2 — \ ^ / 2 \ 
X p \ [r V) 


Vr, 


r'^Tx 


■ 

/ 2 - \ ^ 

— Aul 


- X 

rxj 


n . , 4 

- with V = -Vi + 1/2 > 0. 
r 3 


It should be noticed that V can be rewritten as 


^ V /r\ — d® I (r\ -1^® 


X 


X 




(2,13) 


X y t 


A strong solution to problem 02.121) is defined as follows. 

Definition 2.1 v G L°° ([0, T]; .R))) fl L°° ([0, T]; hT^’°°(/)) with 

r(x, f) = ro(x) + / v{x,s)ds for (x, f) G / x [0, T] 


(2,14) 


satisfying the initial condition O2.I2P 3 is called a strong solution of problem 02.12p in [0,T], 
1) rx{x,t) > 0 for {x,t) G / x [0,T]; 

«;p>/2„eC>([0,r];L7/)); 

j;r€L~(10,r|iffL(|0,fi))) and{p'-y^rfx),, p-'''"V) 6 L“(| 0 ,T];L 2 (/)); 

4) v{0,t) = 0 holds in the sense of W^’°°-trace for t G [0,T]; 

5) 02.12p ^ holds for {x,t) G/ x [0,T], a.e.. 
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We are ready to state the main theorem of this paper. Denote 

- 1, Vx){-,t)\\l^ + (•,t)l |^2 + ||p^/^nt(-,t)||^ 2 . (2.15) 

Theorem 2.2 Let 7 G (4/3, 2), 0 < 0 < 7 / 2 , and p he the Lane-Emden solution satisfying 
CT-(10i). Assume that the compatibility condition n(0, 0) = 0 holds and the initial density 
Po satisfies fll.4p and fl2.9p . There exists a constant S > 0 such that if 

^( 0 ) < 5, 

then the problem fl2.12p admits a unique strong solution in I x [0, cxo) with 

S{t) < t > 0, 

for some constant C independent oft. 

It should be noted that ||p^/^Ut(-, 0) ||^2 is given in terms of the initial data (ro, uq) by equation 

(EHi. 

For any t > 0, since rx{x,t) > 0 for a; G /, r{x,t) dehnes a diffeomorphism from the 
reference domain I to the changing domain {0 < r < R(t)} with the boundary 

R{t) = r (^R, f) . (2-16) 

It also induces a diffeomorphism from the initial domain, to the evolving domain, 

Bn(t){ 0 ), for alH > 0 : 

X ^ 0 G Br^O) r (r^W^Dfi) ^ G BnpfiO), 

where rfi^ is the inverse map of rg dehned in fl2.8p . Here 

Bro{0) ■= {x G : |x| < i?o} and Hj?(t)(0) := {x G : |x| < R{t)}. 

Denote the inverse of the map r{x, t) by TZt for t > 0 so that 

if r = r{x,t) for 0 < r < i?(t), then x = TZt{r). 

For the strong solution (r, u) obtained in Theorem 12.21 we set for 0 < r < R{t) and t > 0, 

p(r,t) = - and u{r,t) = v{x,t) with x = Tlt{r). (2.17) 

r‘‘{x,t)r^[x,t) 

Then the triple {p{r,t),u{r,t), R(t)) {t > 0) dehnes a global strong solution to the free 
boundary problem fll.3p . Furthermore, the strong nonlinear asymptotic stability of the 
Lane-Emden solution can be stated as follows. 

For any l G (0, (27 — 2 — 6 >)/ 8 ], we set 

a = min {7 — 1 + 9, 2(7 — 1 )} — l, (2-18) 

fi = l + (a-t)/(^-e), (2.19) 

<^ = min{l, 2“^/9-I-2“^(/3 — 1) min {1, (7 — 6 ')/a}}. (2.20) 

It should be noted that 1 < /3 < 3 and 0 < a — 9 < (7 — 1) for 4/3 < 7 < 2 and 0 < 6 ^ < 7 / 2 . 












Theorem 2.3 Under the assumptions in Theorem 12.21 the triple {p,u,R{t)) defined by 
fl2.16p and fl2.17p is the unique global strong solution to the free boundary problem fll.ll) 
satisfying R G +oo)). Moreover, the solution satisfies the following estimates: for 

any l G (0, (27 — 2 — 6 *)/ 8 ], I G (0,1) and 6 G [0, 2 — 7 ], there exists positive constants (7^ 


and independent of x and t such that for all t > 0, 

+ ^|r(x,t) — xp < Ct6?(0), X G/, (2-21) 

r 3 max{ 0 , 87 —5+2b} o 1 1 7 ^ 

+ 2 ( 7 -s) ’ j p~^[x)\p{r{x,t),t) — p{x)\^ < C,,(S’{0), X G/, ( 2 . 22 ) 

(1 + (|rj;(x,t) - 1 |^ + |x"V(x,t) - 1 |^) < xg[ 0 ,/], (2.23) 

il + tf/^\u{r{x,t),t)\^ < am, ^ e /, (2.24) 

{1+ ty {\ur{r,t)\‘^ + \r~^u{r,t)\‘^) < C,(S’{0), x G/. (2.25) 

Here a > 0 is given by 


a = min ^ fi — 1, -—-— fi, ^ ^ ~ niax{ 0 , 40 — 4(7 — 1 ) — a}, 


7 + 0 — 0 


4a 


P 




max{0, 40 — 7 — 1} > . 


Remark 2.4 The uniform convergence ofr{x,t) to x is given in fl2.2ip . This implies that 
every spherical surface moving with the fluid converges to the sphere enclosing the same mass 
inside the domain of the Lane-Emden solution with a uniform convergence rate, and the large 
time asymptotic states for the vacuum free boundary problem fll.3l) are determined by the 
initial mass distribution and the total mass. In particular, this also implies the convergence 
of the vacuum boundary: 


\R{t) -R\< Cfil + t)“=G+“-'’)^V(f(0). 

Moreover, a better decay rate for |r — x| away from the vacuum boundary can be given by 

37-2+2(01-9) a 1 

(1 + t) 2 (T'+a-e) P 2 |r(x, t) — Xp < (0), X G [0,/], 0 < / < 1, t>0. 

(See (14.11) for the details.) The uniform convergence of the density and the velocity to those 
of the Lane-Emden solution with uniform convergence rates are given by fl2.22p and fl2.24p , 
respectively. The estimate fl2.22l) yields not only the uniform convergence for large time but 
also the behavior of the density near the vacuum boundary since 7 < 2 . 
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3 Proof of Theorem 12.2 


In this section, we derive some priori estimates under the following a priori assumptions. Let 
n be a strong solution to fl2.12l) in the time interval [0,T] with 


r(a:,t) = ro(x) + / v{x,T)dT, (a:,t)e[0, l]x[0,T], 

Jo 

satisfying the following a priori assumption: 

Ir^; — 1| + |r/a: — 1| < Co for (x,f) G / x [0,T], 


(3.1) 


where eo G ( 0 , 1 / 2 ] is a sufficiently small but fixed constant (indeed, eo is required to be less 
than a constant depending only on 7 ); and 


\'Vx\ + \v/x\ < 1 for (x,f) G / X [0,T]. 

In particular, it holds that 

1/2 < rj;,r/x < 3/2 for (x,f) G / x [0,Tj. 


(3.2) 


(3.3) 


(Indeed, the a priori assumption fl3.ip and 03.21) will be verified in Lemmas 13.61 and [3.101 1 
The following Hardy’s inequalities will be used often to derive the a priori estimates, 
whose proof can be found in [ 8 ] or [T5] . 

Lemma 3.1 (Hardy’s inequalities) Let k > 1 be a given real number and g be a function. 


If g satisfies that x^(g^ + gl)dx < 00 , then it holds that 
/•I/2 |•l/2 


I x^ ‘^g‘^dx< / x^{g‘^ + g‘l)dx < 00 . 

Jo Jo 

Similarly, if ~ + 9 l)dx < 00 , then 


(3.4) 


/ (1 - a;) gdx< (l-x) {g + g,f)dx < 00 . 

n /2 Ji /2 


(3.5) 


3.1 Lower-order estimates 

To derive the lower-order weighted energy estimates, one rewrites equation (12.121) ^ as 




2 — \ 7 ' 
Xp^ 

r 2 


X 


- ZJ (^)x = 


2 — \ ^ / 2 \ 
p \ [r^v)a: 


r 2 


r2r 


■ 

/ 2 - \ ® 

— 4 z/i 

{-A 

J X 

f'xj 


(3.6) 


where fl2.3p has been used. Here z/ = 4 z 2 i /3 -I- z/ 2 . 

We outline the analysis for the lower-order estimates here and give some motivations for 
the proofs. In Lemma 13.21 we derive the decay estimates for the zeroth-order energy: 


x^ < -I- 


r \ 2 

-- 1 ) + r,, 

X / 


ir 


(x, t)dx 
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and the boundedness of the weighted energy: 


\x 


{x, t)dx, 


by constructing suitable nonlinear weighted functionals. In Lemma 13.31 the decay estimates 
for weighted norms of the higher time-derivatives (than those in Lemma 13.2p are derived. 
Further regularity near the vacuum boundary and decay estimates are obtained in Lemmas 
13.41 and 13.51 by showing the boundedness of 




r \ 2 

- - 1) +{r^ 

X / 


ir 


(x, t)dx, 


and the decay of the following weighted quantities: 




r \ 2 

- - 1) + (r^ 

X / 


1)^ 


{x, t)dx, 


x"^ <p v^) -h 


r \ 2 

-- 1 ) +ir. 

X / 


1)^ 


(x, t)dx, 


+ x'^p^vl) {x,t)dx and / p® °‘^'^(x‘^vl + v'^){x,t)dx, 


where a is given by (I2.18p . Those two lemmas play a crucial role to the global regularity 
uniformly up to the vacuum boundary. In the proof of these two lemmas, we use the 
multipliers: 

p"“(l/) (r^ - dp and [ p-^{y){r\)ydy. 


JO JO 

In Lemma 13.51 we also obtain the decay estimates of 


{r{x,t) — xYdx and / x^{rx{x,t) — l^dx. 


With the above lower-order estimates in hand, we can bound 


sup (x^|r3;(x, t) — Ip) and 

X&I 


sup {x^\v^{x,t)\‘^) , 

xGl 


and estimate the decay of 

sup (x|r(x, f) — xp) and sup (x|n(x, f)p) 

xGl xGl 


in Lemma ESI which in particular give the uniform estimates away from the origin. The key 
idea is to integrate fl3.6p both in x and t to derive an ODE for the quantity Z dehned in 
fl3.8ip whose leading term is r^. — 1 to obtain the estimate for x^|r3.(x, t) — Ip. In turn, the 
leading term for dtZ is so that the bound for dtZ yields the bound for x^|na;(x, f)p. In the 
proof of this lemma, the balance of the pressure and gravitational hied plays an important 
role. 
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It worths pointing out the difficulties in the lower-order estimates for the case of density 
dependent viscosities. For example, in the proof of Lemma [3.21 it is shown that the quantity 


{l + t) / < pv"^ + pi 


+ 




(x, t)dx 


<C ^initial data 



2 -7 

X 


(I-!)%(.-It 


dxdsj . 


(3.7) 


In order to obtain the decay of the zeroth-order energy, we show that the double integral on 
the second line of (13.71) can be bounded by the initial data and the term 


C 



p^x 


1 | + 


r 

- - 1 
X 


v\ + p^ \r — x\ \xVx\ 


dxds, 


(3,8) 


using the multiplier — x^ motivated by the virial equations adopted in the study of stellar 
dynamics and equilibriums (cf. [201139]) to detect the detailed balance between the pressure 
and the self-gravitation. The term (13.Sp appears due to the dependence of the viscosities 
on the density. This is in sharp contrast to the case of constant viscosities for which the 
multiplier — x^ matches the viscosities well and the double integral on the second line of 
(13.7p is bounded directly by the initial data. It is quite subtle to bound (13.81) due to the 
degeneracy of the viscosities at vacuum. We choose a cut-off function deliberately whose 
effective length is a small positive number 6 to localize both near the vacuum boundary and 
the origin, so that (13.8p can be bounded by 





x^/F 


(I-i)T(7-U 


dxds 


and other terms. The desired estimates are then obtained by choosing 6 small. In the 
estimates of x^\r^{x,t) — Ip and x^\vx{x,t)\‘^ in Lemma 13.61 we derive an ODE for the 
quantity Z dehned in (I3.8ip . while in the case of constant viscosities, an ODE is also derived 
for the quantity {p{x))~^p{r{x,t),t) which is simpler than that for Z and can be solved 
explicitly in some sense. The ODE for Z is more involved and harder to solve. We have to 
identify the leading terms and estimate the error terms. 

In the following, we give the details and analysis outlined above. 


Lemma 3.2 Let 9 G (0,1]. Suppose that (13.ip holds for suitably small constant Cq. Then, 


SM + 



2 -7 
X 


(^-i)0(7-y 


dxds 


+ / (1 + s) p^ + v'^) dxds < CDi^O), te[0,T], 


(3.9) 


where 


2)(f) =(1 + t) / x^ < pv"^ + p'^ 


(^i)T(.-y 


(x, t)dx 


xV (r.-l)^+(L_i 


(x, t)dx. 
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Proof. The proof consists three steps. 

Step 1. In this step, we prove that for 0 < 0 < 1 and 7 > 4/3, 





(x, t)dx 


<C I ri{x,0)dx + C x^p'^ 





dxds, 


where 


r]{x, t) = -x^pv"^ + 


1 /x\27-2 f 1 

7 — 1 Vr / 



7-1 


+ (iVr 

\ r / * r 7 — 1 


It follows from fl3.6p and the boundary condition fl2.ip that, for any £ > 0 
d 


— <1 (1 + / V{x, t)dx \ + {l + tY I 




— 4z/i (rx^) 


(3.10) 


r^Tx 


= i{l + tY / p{x,t)dx. 


dx 

(3.11) 


Each term in the equation above can be estimated as follows. First, using the Taylor expan¬ 
sion, one may verify that for 7 > 4/3, 


V{x, t) > ^x^px^ ^ xV^ 2 - 1 ) {tx - ly 

2 4 Vx / 


1 


p{x,t) < -x^px^ -I- CYl)x^p'^ 


-- 1 ) +{rx-lf 

X / 


(3.12) 

(3.13) 


provided fl3.ip holds for a suitably small constant eo, where Cf^) is a positive constant 
depending on 7 . Also, 


(r^x) 


- {rv^)x > 3 ^ (yjl + > 


where a = min{2z/i/3, 102 }. Integrate (13.lip over [0,t] to obtain, by virtue of (I3.12p and 
fl3.13p . that for 0 < 6 ^ < 1, 



(^i)T(.-7 


(x, t)dx 


(3.14) 



<C I p(x, 0)dx + Ci I (1 + s)^ W x^ < p®x^ -I- p^ 




dxds. 
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Setting £ = 0 in fl3.14p . leads to 


I pv^ + 


X 


- 1 +{r^-iy 


{x, t)dx 


+ / p^ dxds <C p{x, 0)dx. 


(3.15) 


Letting £ = 1 in fl3.14p and using fl3.15p prove (13.101) . 
Step 2. In this step, we prove that for 9 E (0,1], 


xY 


\x 


dx + 



x'^fp 


< + C 



p^x ( |r^ - 


2(^-l)V(r,-l)' 


- 1 


dxds 


(3.16) 


v\ + |r — x\ \xv^ 


It is easy to check, in view of (13.6p and (12.11) . that 



P 


-1/ / p 


r 4 
X 




Y-x^) 




dxds. 


r X 

X 



x'^ Y Yv) 


It yields from simple calculations that 


(r^ — dx + 4z/i / p' 


-I 



X 


2 \ e 




-(r — X 
r 


X 


=3 


-2 


r 4 

X 


4 (’■“--“) 


X 


r 

2 \ 7 


X 


2 \ 7 


X 


J X 

2 \ 7-1 


r2^ 


Y — Ys 


'xj 1 7 -^ 7 - 

which, together with Taylor’s expansion, gives 


-)%x + 4(-)%,-7(-)V4-, 
r / \r / \r / r 


dx 


dx. 


(3.17) 



P’ 




^3(37-4) f 
>-- X P 


(^i)L(.-7 


dx, 


provided (13.ip holds for small Cq. Simple calculations show that for 9 E (0,1], 


z/ / p 





X 


2 \ ® (^2 


r'^v) 


/ r^r. 


(r^ — x^)^ dx = f x‘^p^po{x, t)dx, 


(3.18) 


(3.19) 


where 


^ 1 / 2 \ ^ 1 

1 / X \ 1 


Po{x,t) = < 


6* \ 


d — 1 \ r^r 


X 


2 \ ^-1 


+ 


d(d-l) 


9 / 2 

x^ , f x^ 
- In 




r^r. 


for 0 < d < 1, 


for d = 1. 
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Moreover, 


T oc^ \ d f 

x'^pvt I- 5 - = — / x^pv 

X J at \x r 


^ ^^dx- x^pv"^ (l + 2 ^\dx 


and 


4z/i / p' 



X 


2 \ e 


/ Lr 


-(r^ — x^) 


dx 


<C / p^x{\r^-l\ + 


- 1 




)Hd. + c//|r-x||x..|rf.. 


It follows from fl3.17p - fl3.2ip that 


^ / <j 3z/xV^ho + ( 3 - ^ ) 




+ ^( 37 - 4 ) /xV 




dx 


<(7 / jp^x ^Ir^; — 1| 4- — — 1 j |x| -I- p® |r — x| |xxa;| -|- x^px^j dx. 


Using the Taylor expansion, shows that for small eo in fl3.ip . 


>- 


2 V r^r 
1 r 


j. 2 r 


r, - 1) + 2 - - 1 

X 


-Uec 


4 \ r'^r.T 


(r. - 1)2 + (- - 1 

Vx 


Note that 


x2p^ (r,^ - 1) + 2 - 1 

Vx 


dx 


= 


(r, - 1)'+ 4 - 1 )' 


dx — 2 



(r — x)^ dx 


= 


(r^ - 1) V 2 - 1 ) dx + 2^ J - 1 ) 

where fl2.3l) has been nsed. We then have, noting fl2.3p again, that 

C~^ J ~ 1 )^ ~ — J x'^P^Vodx 

<cf x^/ (r^ - 1)2 + - ij 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


dx. 


Therefore, (13.161) follows from fl3.22p . the Cauchy inequality and fl3.15p . 
Step 3. In this step, we prove that 


^{t)<C'D{0) + J (1 + s) T-e-lT-n/a y 2:2^61 I 1 ^^ _ ]^|2 _|_ 


- - 1 

X 


dxds, (3.24) 
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where 





^-+ {tx — 1 )^ dxds 

+ J {1 + s) J {x^vl + dxds. 

If (13.241) is true, then (13.91) follows from Grownwall’s inequality. 

It follows from (I3.10j) + fc (l3.16j) with a suitably large constant k that 


^{t) < G 2 )( 0 ) + G 



p^x ( \rx - 1 | + 


- - 1 

X 


v\ + p^ \r — x\ \xVx 


dxds. (3.25) 


Next, we estimate the last term on the right-hand side of fl3.25p . When A > 0, it follows 
from Young’s inequality that for 0 < u; < (7 - 0)/(7 - 1 ), 


(1 + = <! (1 + s)-^ {p^A) 


{p^A) T-® 


7 — 0 


< (1 -|- s) 7-e-‘*'(7-i) A + p'^A. 
Then, for any 0 < a; < (7 — 0)/(7 — 1), 


(1 + s)-^ / 


rx - If + 


- 1 


X 


i-e 


< J (1 -|- S) 7-S-‘^(7-l) I X^P^ 


ir + 


- - 1 

X 


dxds 


dxds 


(3.26) 



X p^ [ \rx - 1 \ + 


- 1 


X 


dxds. 


Let X be a smooth cutoff function satisfying 

X = 0 on [0,1 —2(5], X = 1 on and 0<x<l on [0,1] (3.27) 

for a small constant 6 G (0,1/4] to be determined later. It follows from the Cauchy-Schwarz 
inequality that 

r / r* \ 1 

dx 


(3.28) 


fxp’ 

X (\rx-l\ + 

r 

- - 1 

) n + r — X \xvx\ 


L V 

X 

) J 


<(5^/^(l-I-t) J 

-I- + t)~^ J — lp+ — — 1 dx 

+ (5^^^(1 + ^) y xP^^'^'^ldx + + t)~^ J xp^ k “ ^1^ dx. 

In view of the Hardy inequality fl3.5p and fl2.5p . one gets that 

f XP^~^'^~^'^v^dx < [ ff~^'^~^^v^dx < C f + vl) dx 


' 1/2 


' 1/2 
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and 


XP^ \r — xf dx < [ ^^p^ ^\r — xfdx<C6 f p^ \r — xf dx 

Jl-25 Jl/2 

<C6 [ (^\r — x\‘^ + {vx — 1)'^) dx, 

il/2 


which, together with fl3.28p . imply that 


XP 


x{\rx-l\ + 


- - 1 
X 


m + r — a; \xvx 


dx 


<C5^^^(l+t) j p^ {y‘^ + x^v‘^dx 
+ C5^/\l + t)-^ j (^\rx - Ip + 

Integrating as using 03.261) with u = 1/2, give 


- 1 


X 


dx. 



XP^ 


x{\rx-l\ + 


- - 1 

X 


n + r — x IxVx 


dxds 


<C6^^‘^ / (1 + s) p^ + x'^vl) dxds 


'y — 6 


+ J (1 + s) t-s-(7-i)/2 I ^ \rx — l\‘^ + 

ft r / „ 2’ 


- - 1 

X 


dxds 



+ / / xV ( |rx - Ip + 


- 1 


X 


dxds. 


Using the Cauchy-Schwarz inequality and 02.5p again, one can obtain 



(I-X)p X |r^-l| + 


- - 1 

X 


e-t 1 

<C 6 ^~^ 



'0 


n| + |r — x| Ixn^ 

p® {x^vl + dxds 

e-i)U(.-ir 


dxds 


pt pi— 8 

+ 5x^i+\ f f 

Jo Jo 

<C5X^~^ J j p^ dxds 

pt rl—8 


dxds 


+ C 6 -^ 



2 -7 

X p' 


0 Jo 




dxds 


<U5^"32)(0) + U(53 



2 -7 

X p’ 


e-i)U(.-ir 


dxds, 


(3.29) 


(3.30) 


(3.31) 


where (I3.15p has been used in the last inequality. So, the proof of 03.24p is completed by 
choosing 6 suitably small, and combining (I3.25p . (I3.30p and (13.3ip together. 
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□ 


The following Lemma gives the decay estimates of weighted norms for the higher time- 
derivatives (than those in Lemma [3.2p . 


Lemma 3.3 Let 6 G (0,1]. Suppose that fl3.ip and fl3.2p hold. Then for 0 <t <T, 

{l + t) j [x'^pvf + ^ [v'^ + } {x, t)dx + J {1 + s) j + vl) dxds 

(72)(0) + ^ y Q)dx. 


(3.32) 


< 


Moreover, for t G [0,T], it holds that 


2)i(t) -|- 



2 -7 

X p^ 


(^-i)L(7-y 


dxds 


+ / (1 + s) p^ {x'^vl + v'^ + + '^s) < (72i(0), 


(3.33) 


where 


2i(f) =y 

+P^ 


(rx - 1) 

2 


2 f f' 

2 . I --1 

X 


{x, t)dx + (1 f) / [p{v‘^ + vl) 




(x, t)dx. 


Proof It yields from f fr^ fl3.6p i^ vMx that 


^ I P>ix,t)dx+ I { u 




r. 


2 — \ S /■ 2 \ 
X p \ [r v) 


r^r.. 



r^r.r 


4z/i 


J t 

2 ^ \ 9' 
X p 


2 - \ 9 

X p 


(r Vt)a, - 4z/i — — [{rv)tVt]^ ) dx 


r^ r. 




(rvvt),^ dx 


J t- 


1 L 

2 [r 



X/ t 



r 2 J. 


(27 - l)r^x^ 4- 2(7 - l)rvv:, + ]r—vl 

2r3. 


- / 


X/ J t L 

4^ — S—rx I -I- 2 


dx 

dx 




t 



VVx 


dx =: RHS, 


(3.34) 


where 


$(x, t) 


^ 2 - 2 , 
= 2 ^ P'^t + 




(27 - l)r^x^ -4 2(7 - l)rvvx + '^—vl 
4- ^^vvx . 
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Similar to fl3.12p and fl3.13p . one has 


+ — 


^ 2 -7 

-X p' 


X 




1 


*h(x,t) < -x^pv'l + C(7)x^p’' 


G) 


(3.35) 


(3.36) 


Using the a priori assumptions fl3.ip and fl3.2p . and Cauchy-Schwarz’s inequality, one can 
obtain that 






X p \ [r v) 




{r^vt)x - 4z/i ( [irv)tvt]^ } dx 


>2a / p® {x%l^ + 2v^) dx-C p^ [cq {x^lt + 2^^) + a ^ {x^v^ + v^)] dx 


and 


RHS < ^ J p^ + "^t) ^ J 

where a = min{2z/i/3, 1 ^ 2 }- Here one has used the fact that 0 < 6 ^ < 1 < 7 . Therefore, it 
follows from fl3.34p that 


A 

dt 


$(a:, t)dx + p^ + ’^t) dx < C p^ {x'^xl + dx, 


(3.37) 


provided that eo is small. This, together with fl3.9p . implies 

f ^{x,t)dx + f [ ff [x'^vl^ + v^) {x, s)dxds < ( <h(a:, 0)da: + CS)(0), (3.38) 


which further implies (I3.32p . by using (I3.37p . (I3.9p . (I3.35p . (I3.36p and the fact 0 < 0 < 1 < 7 . 
Finally, fl3.33p is a consequence of fl3.9p and fl3.32p . 

□ 

In the following lemma, we show the improved regularity near the vacuum boundary and 
decay. 


Lemma 3.4 Suppose that (13. ip and (13.2p hold. Let 6 G (0,1], and a and [3 he given respec¬ 
tively in (I2.18P and (I2.19p . If i & (0, (7 — l)/4], then for 0 < f < T, 


D2(t) + j j 
+ j\l + sf-^ j xV 


(^-i)U(.-ir 
(^-i)U(.-7 


dxds 


dxds 


(3.39) 


+ / (1 + s)^ / p^ [x‘^vl + v‘^) (x, s)dxds < C, (S)i(O) + ||rox - lll^cc) , 
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where 


2)2 (t) = / xy' 


2 ^-OL 




(x, t)dx 


+ (1 + / x^f 


(^l)% (.-!)- 


(x, t)dx 


+ (1 + ^ + P 




(x, t)dx. 


Here is a constant continuously depending on l, but not on t. 


Proof. The proof consists of three steps. Recall that a and ft are given respectively in (12.181) 
and fl2.19j) . that is, 

a = niin {7 — 1 + 6, 2(7 — 1 )} — t and ft = 1 + {a — L)/{'y — 6). 

Step 1. In this step, we prove that for any cu > 0 and k > 1, 


^ 2 - 0 -a 


(i-i)T(.-if 


(x, t)dx 



+ I I xV“ 


X 


--1 +(r,-l)^ 


dxds 


<C ( 2 ) 1 ( 0 ) + llrox - l||ioo) +Cu: / (1 + s)'" / + x^vl)dxds 


>0 


Cuj W (1 + s) 

Jo 


-K / ^2-e-a 


r \ 2 

-- 1 ) +{+ 
X / 


1)^ 


dxds. 


(3.40) 


It should be noted that k > 1, which ensures the last line in fl3.40p behaves well according 
to the Gronwall inequality. 

Multiplying (13.Op by P~‘^(y) — y^)ydy and integrating the resulting equation with 

respect to the spatial variable, one has, with the help of the integration by parts and boundary 
condition ( 12 .ip . that 


P’ 


px 


P "(y) -y^)dy 




2 \ ^ 


(r^v) 


p{^) p "(y){^"-y")ydydx-u J P^{r^-x^)jx 


ry^^ry^ j ry^Ztryi 


+ 4z/i / p® 


X 


2 \ 0 r 


r^r 


V 

r JO 


P “( 2 /) -y^)dy 


dx. 


As the derivation of (I3.16p . one can get 


x2/-“ 


(^-i)A(.-y 


dx + 



x2p^-“ 


(^i)%(.-7 


< 


X 
3 


dxds 


C f x^p^ “ — 1 ') + (r^; — 1)^ {x,0)dx + C'y2 [ \J^i\ds, 

J l\x J J io 


(3.41) 
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where 


Li = 


_x 




L 2 = 



p -y^)ydyj dx, 


L 3 =4z/i / p' 


—I 



X L 

2 ^ ^ r 


dx, 


r^r 


P “( 2 /) {r^ -y^)dy 


dx. 


In the above integration by parts, we should notice that O<a<0 + (7 — 1) which ensures 
that p^~°‘ is integrable on [0,1] due to 02.51) . so that the above integrations by parts are 
legitimate. 

Next, we estimate Li, L 2 and L 3 . For Li, it follows from the Cauchy inequality that for 
any cu > 0 , 


|.hi| < C / pvt 


P “ 1 / (k/y - 1| + Vy - l\)dy 


dx 


<Cu-^ / p^-^^-^'^vldx + UJ 


Due to 02.51) and 03.5p . one has 


P “y (k/y- 1| + Vy - l\)dy 


dx. (3.42) 


r /• 1/2 /•! 

/ <C ^v^dx + C + vl)dx 

J Jo Jl/2 

^1/2 p1 j‘ 

<C / p^v^dx + C / p^{vl + x^vl^jJjdx < C p^{v^ + x‘^Vf^)dx. 

Jo Jl /2 J 


(3.43) 


Due to Holder’s inequality, 0 < 1 and a < 2(7 — 1) (which implies a + 9 < 27 — 1), and 
02 .5p . one obtains 


pi-e+i 


P "d i\r/y-M + \ry - l\)dy 


dx 



< / pi e+7 I / p 7 ( / pT' »y‘^(^\r/y-l\ + \ry-l\)dy]dx 


<C / V(k/ 2 / - 1| + K - l\)dy. 

Jo 

It yields from 03.421) - 03.44p that for any a; > 0, 

[ |Fi|ds [ [ ff {Vg + x‘^v'^,,^)dxds 


(3.44) 


+ Cu 



x^p^-^ 


(I-i)N(.-7 


dxds. 


(3.46) 
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Rewrite L 2 as 


Lo — 



(p -y'^)dydx 


/ -p-[r^-y\dy 


dx —: L 21 + Pv 22 - 


It follows from that fl2.3p . and the Cauchy and Holder inequalities that 
/• 1/2 


|^2l| - ^ j ^ ^ (1^^ - 1| + 

»l/2 

_ 1 
X 


- - 1 

X 


<CJ^ + 

<C j x^ - 1|^ + 
and for any ca > 0 , 


y^\r - y\dy 
1/2 


dx 



- 1 


X 


dx + C I x~^ 

Jo 

/•1/2 

dx<C x^p^ ( Ira; — 1 |^ + 


y^dy] ( / \r-y\^dy]dx (3.46) 


- 1 


X 


dx, 


|i:22|<C'/ /C-“(|ra;-l| + 

7i/2 


Jl/2 

<u: f 

Jl/2 


- - 1 

X 


j |r — x|da: 



p-V k,-i| + 


- -1 
y 



dy ) dx 


' 1/2 


/C "(r — x^dx 


(3.47) 


+ Cuj-^ / p^+“ ( / p-^-^‘^y‘^dy 

Jl/2 



p^y |^|ry-l| + 

By virtue of (12.51) . (13.51) and a < 2(7 — 1), one has 
/•i /-i 

pl-^{j — x^dx <C / p'^“"+2(7-1) _ 3,^2 _ ]^^ 2 j 


-1 


y 



dy dx. 


' 1/2 


' 1/2 

»i 


<C / p'^ [{r — x)"^ + x'^{rx — l)"^] 

Jl/2 


and 


' 1/2 


^ 7 +« / p-l-2ay2^y \dx<C. 


It then yields from (I3.46P and (I3.47P that for any cn > 0, 



(^-i)H(.-ir 
(^-i)C,.-if 


dxds 


(3.48) 


dxds. 
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For L 3 , it holds that 


iLsl <C J ^ \{r^ -x^),\dx + C J ff 
<C J {x\rx — 1\ + \r — x\) dx 
+ C [ p^x~‘^{x\v^\ + \v\) [ p~'^{r^-y^)ydy 


p ^{r -y )ydy 


dx 


(3.49) 


dx —: L 31 + L 32 . 


It follows from (13.4^ . (13.dh . (I2.5p . 6* — a > 1 — 7 and a < 2(7 — 1) that 

r /•I/2 rl rl /2 

/ p^~°‘v‘^dx<C / v^dx+ / p^~'^v^dx<C / x‘^{v^ + vl)dx 

J Jo Jl/2 Jo 


-1 


(3.50) 


+(7 / /-“+2(7-1)(^2 ^ ^ 3.2^2)^^^ 

Jl/2 J 

which, together with the Cauchy inequality, gives that for any cn > 0, 

L31 <Ca;(l + t)'^ J + x^v 1 )dx 

+ + f)-" J xV"" - 1 ) ' + (r^ - 1)' 

Clearly, L 32 can be bounded by 


(3.51) 


dx. 


1 . 1/2 

L 32 <c / x~‘^{x\vx\ + |n|) 

Jo 


p-^{r^-x^)- / [p--)y[r^ - y^)dy 
Jo 


[ p^{x\vY + 1^1) 

f p V fkj/-1| + 

r 

- - 1 

I dy 

/ 1/2 

Jo V 

y 

J 


dx 

dx =: L 321 + F 322 . 


For L 321 , it follows from fl2.3p and the Cauchy and Holder inequalities that 
1.1/2 1.1/2 


F 321 < C / {x\vx\ + \v\)\r — x\dx Jr C / x~^{x\vx\ + \v\) 

Jo Jo 

- 1/2 1.1/2 


v V -y\dy 


dx 


^1/2 pi 

<C / [p^{x‘^vl + + p'^{r — xY] dx + C / ( / y^dy \ ( / {r — y^dy ] dx 

pl/2 

<C / [/i^(x^n^ + n^) +//^(r — x)^] dx. 

Jo 

For L 322 , it follows from the Cauchy and Holder inequalities that for any cu > 0, 

L322 <1^(1 + ty [ + v^dx 

Jl/2 

+ Cu}-^{1 + t)-^ [ {\ry - 1|2 + \r/y - 1|^) dy. 
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Here one has used the fact that a < 2(7 — 1), which implies 


/ nx 

^0 --6/-a_2, 


J 1/2^ [Jo p y 

So, it holds that for any cu > 0, 

L32 <a;(l + ty J + x^vl)dx + ^ J dx 

+ Cu~^{l + ty^ J ^—— 1 j + (ra; — 1)"* dx. 

This, together with (I3.49p and (I3.5ip . implies 

/•£ pt p 

/ 7 ^ 6>/_2 I _ 2 _ 2 ', 


\L 3 \ds <Cu / (1 + s)'^ / + x^v^dxds 



(3.52) 


—K I _2 —O—ct 




dxds. 


The estimate fl3.40p follows from fl3.45p . fl3.48p . fl3.52p . fl3.33p . and the fact that r(0, t) = 0 
(which implies ||r — xjl^^oo < Wr^ — 

Step 2. In this step, we prove that for any k > 1, 


^r(t) < C'Di(O) + C f (1 + f xV-“ Q-iy + y - 1) 


dxds 


i-e 


+ C / (1 + s) (1 + s) 

Jo 


ft -1 / ^2 -e 


xV |ra,-l|^ + 


- -1 

X 


dxds 


+ C(l + f)^ J x‘^pv^{x,t)dx + C J (1 + s)'' J p^{x‘^v^ + v^dxds, (3.53) 


where 


^^{t) ={l+ty J x"^ l^pv"^ + y {~~y +(ra:-l) 

+ {i + tr-^jxy^ y^-iy + y-i) 

+ f {1 + sy f p^ {x^vl + v^) dxds 


(x, t)dx 


(x, t)dx 


+ / (1 + sy-^ / xV 


(^-i)T(.-7 


dxds. 
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It follows from a suitable combination of (1 + ^ (\?).22\\ ds with fl3.14p ioi i = k that 


^r(t) < C2)i(0) + C / (1 + s)^-2 / 


+ C (1 + s)^ / pxllr^; —1| + 




dxds 


- 1 


X 


|n| + |r — x| \xVx 


dxds 


+ C(l + t)'^ ^ / x^pv‘^{x,t)dx + C I (1 + s)^ ^ / p^v^dxds, 


(3.54) 


where 03.231) and 03.33P have been used. Each term on the right-hand side of 03.54p can be 
estimated as follows. It follows from the Young inequality that for any cu > 0, 


[\i + sy-^ [xV 


Jo J 

\x J 


[\i + sy-^ [ xV 


Jo J 

\x / 


dxds 

dxds 


(3.55) 


[ (1 + 

[ x2/-“ 

(--l)' + (r.-l)^ 

Jo 

/ 

Vx / 


dxds, 


due to 

(1 + s)^~'^ff = ^(1 -|- 5) l-e + a + a _|_ s) 7-e + a 

Let X be a smooth cutoff function satisfying 03.27p . Similar to fl3.29p and 03.26p . one obtains 


(1 + 1 ) 


xf 


x[\r^-l\ + 


X 


n +\r — x\ \xvx 


dx 


+ ty j jf {v^ + x\l) dx 

+ + ty-^ j (^|r, - 1|2 


(3.56) 


and 


k- 2 / 2-e+(7-l)/2 


(1 + ty-^ / X^p' 


|r.-lp + 


- - 1 

X 


<(1 + ty ^ 7-e-'(7-i)/2 / x^p^ ( Ir^ — 1|^ + 


■ 7 —6 


2^6» 


- - 1 

X 


dx 

2 


dx 


+ {i + ty-^ / xV( |r.-ir + 


-1 


X 


- - 1 

X 

dx. 


dx 


(3.57) 
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It thus yields from fl3.56p and fl3.57p that 


[\i+sr-^ [ x(P 

“T 

1 

+ 

r 

- - 1 

] X + r — X XXa, 

Jo J 


X 

/ 


dxds 


y (1 + s)^ j {y^ + x^vl) dxds 

+ 06^^"^ J (1 + s)''“^“7-e-(7-i)/2 J x'^p^ — 1\'^ + 

+ ^ (1 + 5)"-' y (^|r, - 1+ 

Similar to fl3.3ip . one can get 


- 1 


- - 1 

X 


X 


dxds. 


(3.58) 


dxds 


(l + ,)-i {l-x)p^ X |r.-l| + 


- - 1 

X 


n +\r — x\ 


dxds 


<(75^-1 ^ I {1 + s)'^ ^ I p^ dxds 


(3.59) 


+ C6-^ / (1 + s)"-^ / x^fP 


dxds. 


Now, the proof of (I3.53p is completed by choosing u and <5 suitably small, and combining 
(ESP, (13351) . (13581) and (13591) . 

Step 3. In this step, we prove that 


^i{t) <C (2)^(0)+ ||ro.-l||i.), 
where (3 is defined by (I2.19p . and 


(3.60) 


2-6/-a 


^2{t) =^l{t) + / X^P 



(^-i)%(.-i) 

( 3 i)N(.-i) 


(x, t)dx 


dxds for any k > 1. 


It follows from (I3.40p and (I3.53P that for any k > 1, 

^ (®l(0) + Ikox — l|lioo) + C f (1 + s) ^ 


+ (1 + s) '"j j xV “ ~ 


dxds 


'y — 9 


+ C / (1 + s) 7-e-(7-i)/2(l + s) 

Jo 


K-1 / „2-6» 


2 -y / I 1 12 I 

X p |r,j, - 1| + 


X 


dxds 


+ C(l + t)^ J x^pv^{x,t)dx + C J (1 + s)'^ J p^{x^v^ + vydxds. (3.61) 
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When /9 < 2, fl3.60p follows from follows from fl3.6ip with k = (], fl3.33p . 9 < 1, and the 
Grownwall inequality. When (3 > 2, we have 

l-a/i'y-e) < -9). 

Then, it follows from (13.6ip with k = 2, (I3.33p . 9 < 1, and the Grownwall inequality that 
^l(t) <G (Di(0) + ||ro. - . (3.62) 


Due to 6* < 1, one has 


(3 = 13 -^ < IH-r < 3. 


7 — 0 ' 7 — 1 

So, 03.601) follows from 03.6ip with k = (3, 03.621) . 0 < 1, and the Grownwall inequality. 
As an immediate consequence of 03.39p . we have the following Lemma. 


□ 


Lemma 3.5 Suppose that 03.ip and 03.21) hold. Let 9 G (0,7/2], and a, {3 and be given 
respectively in 02.181) . 02.191) and 02.201) . If l ^ (0, (27 — 2 — 0)/8], then for 0 < t <T, 


2^-2+a-9 


(l + t) T+Q-e ^ J (r{x,t) — xfdx + {1 + t)'y+<^-o^ J x^{rx{x,t) — lydx 
+ (1 + ty j {x‘^pv‘1 + p^^x^vl + (x, t)dx 

+ [ y + + vl)dxds < C, (S)i(O) + llroa; - lH^oc) , 


(3.63) 


{1 + ty J (v^ + x'^p^vl) {x,t)dx + {l + t)^2' J °‘^'^{x‘^vl + vy{x,t)dx 

< Cl (2)i(0) + llrox — l||^oo + llw(') 0)||^oo) . 

Proof. The proof consists of three steps. 

Step 1. In this step, we prove that 

(l + t) 9 j (r(x,t) - xfdx < Cl (2)i(0) + \\ro^ - 1||J<^) , 

(1 + J x'^yyx, t) - ifdx < Cl (2?i( 0) + ||roa; - l||5^oc) , 

(1 + tY J {x^pvl + p^{x^vl + vY) (x, t)dx 

+ {l + sf J p\x^vl^ + vl)dxds < Cl (S)i(O) + ||rox - lH^oc) . 


(3.64) 


(3.65) 

(3.66) 


(3.67) 
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It follows from Hardy’s inequalities fl3.4p and fl3.5p . fl2.5p . 9 < a and the Holder inequality 
that 


<C x^p^~°‘ ((r — xY + {r^ — 1)^) 

((r — xY + Yx — 1)^) dx^ 


2-1 

7 + q : —6 


27 —2+ q ; —0 
7 + q : —6 


X / X 


Thus, 


Clearly, 


X {Tx — 1) dx <C / X p “(fa; — 1) dx 


7 

7+0: —0 


x^p'^Y'x — l)^dx 


0 — 0 
7 + 0-0 


With fl3.39p . one may easily obtain fl3.67p by virtue of fl3.37p . fl3.35p . fl3.36p and fl3.33p . 
Step 2. In this step, we prove that 


(1 + tY J [{x^Yx - 1)^ + {r - x)^)] (x, t)dx 
+ [ {l + sY [ p^~°‘{x^vl + vYdxds < C, (S)i(O) + ||rox - l||5,oc) 


(3.68) 


Multiplying equation (13. 6 p by p °‘{y){r'^v)ydy and integrating the product with respect 
to spatial variable, one obtains, following the derivation of fl3.10p . that 


(1 + ty / x^p^-" 


— 1 j + {tx — 1)^ (x, t)dx + [ {1 + sY f p^ dxds 


< C ||ro. - lllioo +CY^j\l + sYmds + C J xV" (^ - 1 )' + - 1) 


dxds, 

(3.69) 


where 


X 


/ p '^)ydy dx, 


X 


^2 = / ^ p '^r V - p ^{r v)ydy 


K 3 =4z/i 



dx. 


p °'{r^v)ydy — p “r^x 


dx. 


The estimate for Ki can be obtained in a similar way as the derivation of fl3.45l) . fl3.48p and 
fl3.52p . Clearly, 


(1 + s)^|J’7i| <Ccx ^ / {1 + sY p^v^dxds 


(3.70) 


+ Ccx / (1 + s)^ / p® “(x^ + x^x^)dxds 










































for any w > 0 . To deal with i^2, one notes that 


(1 + ty / p'*' " ix\rx — 1| + |r — xl) \v\dx 

Jl/2 

<(1 + tf [ p^-^V^dx + (1 + [ p27-«-e _ ^|2 ^ 1^ _ ^|2^ 

Jl/2 Jl/2 

<C(1 + t)^ [ p^{x‘^vl + V^)dx + (1 + [ p27-a-e (^x‘^\r^ _ 1|2 _|_ 1^ _ 


' 1/2 

This implies that if 6^ + a < 7, 


'1/2 


{1 + ty p^ {x\rx — 1| + |r — x|) \v\dx 

Jl/2 


<C{1 + ty J p^{xvl + v^dx + C{1 + ty ^ J p'^ {x'^\rx — + |r — xp) dx] 


and if 0 + a > 7, 


(1 + ty / p'’' " {x\rx — 1| + |r — x|) \v\dx 

Jl/2 

<C(1 + ty J p^{xvl + v^dx + ^ y P^~'^ — 1|^ + |r — x|^) dx 

+ C(1 + f p^ {x‘^\rx — if + \r — xf) dx, 


due to 


Then, one can obtain 


a 


-y — 9 


{2<y-y)<y-i. 


[ (1 + s)'^|it' 2 |ds <C f (1 + s)^ [ p^{xvl + v^dxds 

Jo Jo J 

+ Cy y p"^”" (x^lr^; — Ip + |r — x|^) dxds 
+ Cy (1 + s)^“^ y fP [x‘^\rx — if + \r — xf) dxds. 

Note that for any w > 0, 

[ {l + sy\Kyds <Cu~^ [ (1 + s)*^ f ff {xvl + v‘^)dxds 

Jo Jo J 

+ 0 ; y {I + sy J p^~°‘{xvl + vydxds. 

So, follows from (EJl, flT89|l and flMjl - flTTD . 


(3.71) 


(3.72) 
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Step 3. In this step, we prove that 

(1 + t)^ j + v‘^){x,t)dx < C, (S)i(O) + ||rox - lH^oo + ||yx(-, 0)||Joo) , (3.73) 

{1 + tf J + x^p%l){x,t)dx < C, (S)i(O) + ||rox - lH^oo + ||wx(-, 0)||^oo) • (3.74) 

Notice that 


(l + t)"^ I p^-^/\x^vl + v^){x,t)dx 


= f + v'^){x,0)dx + [ (l + s)^ 2 '' ^ [ p^ + v'^)dxds 


'0 


+ 2^ (l + s)^2 J ^^^{x'^VxVsx + vVs)dxds 

<C ||nx(-, 0 ) 11^00 +C y (1 + J + x^vlx + vl)dxds 

+ C J {1 + sY J ff~°‘{x‘^vl + vYdxds. 

Then, fl3.73p follows from fl3.39p . fl3.67p and fl3.68p . Similarly, one can show that 



(1 + tf / p\x^vl + vY{x,t)dx < C, (Di(0) + ||rox - lH^oc + ||'y2;(-, 0)||^„e) • 


This, together with (12.5p . (13.4p . (13.5p and 2(7 — 1) > 7/2 > 9, gives (I3.74p . 

□ 

Suppose that fl3.ll) and (13.21) hold. Let 9 G ( 0 , 7 / 2 ], and a and (3 be given respectively in 
fl2.18p and fl2.19p . Then, for any l G (0, (27 — 2 — 0)/8] and / G (0,1), one has 

Lemma 3.6 Suppose that fl3.ip and (13.21) hold. Let 9 G ( 0 , 7 / 2 ], and a, {3 and ^ be given 
respectively in (12.181) . fl2.19p and (12.201) . If i E (0, (27 —2 — 6 ')/ 8 ], then for {x,t) G Jx [0,T], 

+ x^\rx{x,t) - 1^ + x^\vx{x,t)\^ < C, (S)i(O) + Ijrox - l||^oo + ||wx(-, 0)||^oo) . (3.75) 

Proof. The proof consists of three steps. 

Step 1. In this step, we prove 

(1 + t)\+«-« ^x|r(x, t) - xp < C, (Di(0) + \\rox - l||5,cc) , (3.76) 

(l+t)^"^”“^°’ ^^~‘^^'^~^^~°‘^xv^{x,t) < C, (2)i(0) + |]roa; - l||ioc + ||'yx(-,0)||^oo) . (3.77) 
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Clearly, fl3.76p follows from fl3.63p and 

pX 

x{r{x,t)-xf= {y{r{y,t)-yf)dy 


< J {r — xy{x, t)dx + 2 (r — x^dxj 
Similarly, for x G /, it holds that 


y 

\ 1/2 


\ 1/2 

x‘^{rx — l^dx j , x e /. 


(3.78) 


xv^{x,t) = J {yv‘^)^dy < J v‘^{x,t)dx + 2 ^v'^dx^ °‘^‘^x‘^vldx^ 

It follows from fl2.5p and fl3.5l) that 

/ <C V^dx + C p«/2-0+2(7-l)(y2 ^2^^^ 

J Jo Jl /2 

2a-49+4(7-l) 49-4(7-l)-ct 

<C J v^dx + C (^j x^p^vldx^ x^ff~°'^‘^vldx 

if 6*>7 — l + a/4, and if^^y — l + a/4, 

f < C f v‘^dx + C [ x^p^v^dx. 


Then, fl3.77p follows from fl3.64p . 

Step 2. In this step, we prove 

x^\r^{x,t) - 1|^ < C, (Di( 0) + llrox - l\\l^ + ||wa;(-, 0)||5,oo) • 
Rewrite equation fl3.6p as 


Set 



X 


2 \ 7 


r^r. 


- -I (/!"). 


7-6/ 


9 ■ 

Integrate equation (I3.80p over [x, 1] to get 

^2 \ 7-6 


9 


-2 


dtZ + ^p'^ \x)‘ 


j,2j. 


Z + p^-\x) 


3.2 \ 7-6 ±ae_4 


r'^r^ 


- dtY{x,t) -'^i2i{x,t), 


(3.79) 


(3.80) 


v{y,t)dy. (3.81) 


(3.82) 


i=l 
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where Y and satisfies the following estimates: 


u 




= p \x) I 



dy, 


4:Ul6 


-4 


dy, 


Rewrite Z as 

^ - 1 ) + 

where 

YiI-28 


Zi{x,t) = ^ 




Under the a priori assumptions (13. ip . we can see that the leading term of Z i 
Note that 


\x) 


X 


r2^ 


2 \ T'-® ±a0_4 




where 


\£,{x,t)\<Cx ^{x)\r{x,t) — x\. 

Then, fl3.ip and the Taylor expansion imply that 

3 

dtZ + p^~^{x)a{x, t)Z = —dtY{x, t) — ^ S.i{x, t), 


2 = 1 


where 

7 /( 2 © < a{x,t) < 2'^lv, 

\2^z{x,t)\<Cx~^f~\x)\r{x,t)-x\+C^~^\x) f p{y)\v{y,t)\dy, 

J X 

due to the smallness of eo- Integrate fl3.84p on [0,t] to get 

3 rt 


Z{x, t) < Z{x,0) + C sup \Y{x, s)| + 

sG[0,t] 


^ u,p|_pT-»(,) r 

fci Jo I. J- 


a{x, T)dT 


(3.83) 


u{l - Tx). 


(3.84) 


\2oi{x, s)|(is. 
(3.85) 
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Next, one needs to estimate the terms on the right-hand side of fl3.85p . It follows from 
a > 9, 9 < 7/2 and fl2.5p that 


\Y{x,t)\<Cp (x) y p (y) {y\ry - 1\ + \r - y\) dy 


1/2 / .1 


p^ “(2/) {y^\ry - Ip + |r - 


SCp-^x) ( / v-yr“(y)didj 

\ 1/2 

( / p^-^{y){y‘^\ry-l\^+ \r-y\^)dy\ , 



1/2 


\2.i{x,t)\<Cp \x) / y ^p<{y){y\ry-l\ + \r-y\)dy 

J X 

1/2 / pi 


( [' (sV, - Ip + 


“(2/) (2 /Vy-l|^ + |i’ 

1/2 


\S^ 2 {x,t) \ < Cp'^ \x) 

|'23(ai,t)| < Cp^~\x) 
Notice that 



X y'^pv^dy^ +x ^\\{xv^){-,t)\\L^ 

x~^\r{x,t) — x\ + (J y‘^pv^dy\ 


P^ ^(x) / exp < — ^(x) / a(x,T)dTi ds < C 


and 


(3.86) 



1/2 


(3.87) 


(3.88) 

(3.89) 


\Z{x,t)\ < C\rx{x,t) — l\ +Cx ^\r{x,t) — x\ + Cx (^J V^px^dy^ . (3.90) 

Then, fl3T9]l follows from (ESS]), fiTMIl - fiTMl . (EEH]), fiTTOll and fiTTTD . 

Step 3. In this step, we prove 

x^\v,;{x,t)\‘^ < C, (Di(0) \\rox - llljoc + ||iix(-, 0)115^0,) . (3.91) 

It follows from fl3.83l) . fl3.ll) and (13.21) that 
\vx{x,t)\ < C\dtZ{x,t)\ + C\dtZi{x,t)\ 

and 

\dtZi{x,t)\ <Cx~^\v{x,t)\ + Ceo\vx{x,t)\ + Cp~\x) j p{y){\vt{y,t)\ + \v{y,t)\)dy 

J X 

<Cx~^\v{x,t)\ + CeQ\vx{x,t)\ + Cx~^^‘^ (^j y‘^p{v^ + vl)d'i^ 
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Due to the smallness of eo, one has 


\v-rXx-,t)\ < C\dtZ{x,t)\ + Cx ^\v{x,t)\+Cx ( I y‘^p{v‘^ + v^)dy 


»i 


1/2 


Note that 


(3.92) 


\dtY{x,t)\ <Cp / y P {y){y\vy\ + \v\)dy 
J X 


<Cp 


T.-e 




0 

1/2 


{y){yVy + v)dy 


1/2 


(3.93) 


<Cx \y){y^vl + v‘^)dy^ , 


due to P < 7/2 < 2(7 - 1). Then, fim]) follows from flT^ . fiTMl . flT8711 - fl3:mi|l . fl3T6|l . 
fl3T7)l . fl3T9|l . (ESnD, (1031) and flTMll . 

□ 


3.2 Higher-order estimates 

For the higher-order estimates, we set 
^2 


X 


£1 = -^ -1 


r2r 


(3.94) 


Then for any positive constant /c, 
„2 \ 


P 


X 


r^Vn 


- 1 


pM (1 + £3)" - 1 = kp^ (1 + £3)'^-^ £3.. + ^^(x, t)- (3.95) 


where satishes the following estimate: 

\^k{x,t)\ < Cxp^~£~^^ {\r^ - 1| + \r/x - 1|), 
due to (12.31) . Using above notations and (12.31) . one can rewrite (13.61) as 


(3.96) 


u 


4i/-| g -I 

(^) (p/ (1 + £3)'"-' n. + ^e) + (^) [7P^ (1 + £3)^"' £3. + 


X 


Iri/ 


X\‘ 


r J 


X 


P\z) D + -7 - 1 (p^) 


/T \ “^^“1 


(3.97) 


7 \x 


The principal part on the left-hand side of fl3.97p is 

vp^£l^t + iP^O-x = z^(/£3a:)t + 7p'^“^(/£3a;), 
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which can be understood as a damped transport operator for p^ 0 .x with a degenerate damp¬ 
ing coefficient . This is an interplay between the viscosity and the pressure. This 
structure leads to desirable estimates on the derivatives of H, for example, the bound for 

J + j J p^'^~^~^£lldxds. 

The bounds for the weighted norms of the derivatives of £J in turn yield the bounds for those 
of the derivatives of r and v, which are given by Lemma [3.71 stated below. 

In Lemma 13.81 one can use Lemma 13.71 to bound 


p^^ ^ [{{r/x)xf + rl^] {x,t)dx 


and to estimate the decay of 


+rlx] {x,t)dx, 0 < / < 1 . 


Based on this and other estimates, we can estimate the decay of 


{x,t)dx and / \{,{y/x)x) + {x,t)dx, 0 < / < 1 


Putting those together yields the desired decay estimate for 


{Vx, v/x, Vx - 1, r/x - l)(-,i 




, 0 < / < 1 , 


which, together with the supreme norm estimates for — 1 and Vx away from the origin 
given in Lemmaverifies the a priori assumptions fl3.ip and fl3.2p and closes the bootstrap 
argument. 

We now carry out the strategy outlined above. 


Lemma 3.7 Suppose that fl3.1l) and fl3.2p hold. Let w{x) be a smooth function on [0,1] 
satisfying w{x) > 0 and w\x) < 0 on [0,1], and to(l) = 0. Then, 


J w [4(r/x — 1 )^ -1- (ux — 1)^] dx <2 J w£fdx, 

J w [{v/xY + {vxY] dx < C J w{ 0 .tYdx, 

j w [{{r/x)xf + [vxxf] dx <C j w{£lxfdx, 

J w [{{v/x)xf + [vxxf] dx<C J w dx. 

Proof. It follows from fl3.ll) that 

n = - 2 {r/x - 1 ) - {tx - 1 ) + 0 {l)eo {\r/x - 1 | + \rx - 1 |), 


(3.98) 

(3.99) 

(3.100) 

(3.101) 
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which, together with the integration by parts, 'r(0,t) = 0 and w\x) < 0 , gives 
> / » [4(r/x - 1)^ + (r. - 1)^ d. - 2 /(.-„).(r - 
— Ceo J w \{r/x — 1)^ + {r^ — 1)^] dx (3.102) 

y w [4(r/a: - 1)^ + {r^ - 1)^] dx, 

due to {x~^w)x < 0 and the smallness of Cq. This proves fl3.98p . Similarly, fl3.99p follows 
from 


Oi = -(1 + 0 ) {2v/r + v^/r^). 

For (13.loop , note that 

0 ^ = -(1 + 0 ) [2(r/x)(r/x)^ + r^^/r^] . 

Thus, 

j w(0x)^ >2c j w [2(r/a;) X ^^0 ^ ^ [^((^/ ^^X^ H~ ij"XX^ j 

>c J w [{{r/x),^f + (r 

XX )^] dx, 

due to the smallness of Cq, and 

4 J w{r/x)^r^^dx=% J w {{r/x)^f dx - 2 j {xw)a,{{r/x)^f dx 

=6 j w{{r/x)xY dx —2 J xw^{{r/x)^Y dx >Q. 

Similarly, one can prove fl3.10ip . 


(3.103) 


(3.104) 


□ 


Lemma 3.8 Suppose that fl3.ip and (13.21) hold. Let 6 G (0,7/2], and a and /3 be given 
respectively in (12.181) and (12.191) . If l E (0, (27 — 2 — 6)/S] and I G (0,1), then for t G [0,T], 


j ^ W/x)^f E-rl^] {x,t)dx < C(f(0), 

(1 + tf-^ [ [{{r/x)^f + rl^] (x, f)dx < C^S{Q), 
Jo 

y (1 + s)^“^ y [vl + iv/xY) dxds < C,(S’{0). 
Proof. The proof consists three steps. 


(3.105) 

(3.106) 

(3.107) 
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Step 1. In this step, we prove 

[ p'^'^~^{x)0.l{x,t)dx + f [ p^'^~^~^O.ldxds < CS’{0). 


(3.108) 


This, together with flO.lOOp . gives fl3.105p . 

It follows from fl3.97p . fl2.3p and the Cauchy inequality that 


J p'^'^~^\£l^{x,t)\^dx + J J 10.^1“^dx 

<0(0(0) + ^ y d(^) — 1)^ + (r — (x, t)dx + C 

py+l-e _ 3.^2 ^2j 



p^-^^-^^v^dxds (3.109) 



tc 

Notice that 

f p^~^'^~^'>v'^dx <C j x'^{v'^ + vl) + C [ + vl)dx 

J Jo Jl/2 

<C J p^{v'^ + x‘^vl)dx, 

where (13.4p . (13.5p and (12. 5 p have been used. Then. (l3.108p follows from (I3.33p and 6^ < 1 < 7. 

Step 2. In this step, we prove fl3.106p . Let "0 be a smooth cut-off function defined on 
[0,1] satisfying 

0 = 1 on [0,1 — /], 0 = 0 on [1 —//2,1] and — 8 //< 0'(x) < 0 on [0,1] (3.110) 
for any fixed constant / G (0,1). As shown in Step 1, one can obatin 


[ 0p27-l-2^ 


9 dt 


In ( 


dx+ 


^ 37 - 1-0 




<Cy (f{v^ + x^v1 + vl)dx + C j [x^(r3; — 1)^-f (r — x)^] dx, (3.111) 

which, together with fl3.39p . fl3.63p and fl3.108p . implies that 

(1 -h t)^~^ J i>p‘^'^~^{x)£ll{x, t)dx + J (1 + s)^~^ J 'i/jp^'^~^~^O.ldxds < C^(g‘{ 0 ). (3.112) 

(If 0 > 2, we first show the time decay with rate —1. With this, we can then show the time 
decay with rate 1 — 0. Indeed, this technique has been used in Step 3 of the proof of Lemma 

E31) 

Squaring fl3.97p . multiplying the resulting equation by and integrating the 

product over [0,1], one obatins that 

[ {I + s)^~^ [ 'ijjp^"'~^~^£llgdxds < C / (IT s)^“^ [ ip p^'^~^~^O.ldxds 


+ 0(1 + 3} 


3-1 


ipp 


® ^x‘^(rx — 1)^ + (r — x)^ + ^1 + x‘^v‘1 + n^] dxds, (3.113) 
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where fl3.2p has been used. This, together with fl3.112p . fl3.39p . fl3.63p and 0 < 1 < 7, gives 


^ (1 + sf-^ j < a,/^(0). 

So, it follows from fl3.100p fl3.10ip . fl3.112p and fl3.114p that 
(1 + tf-^ J [{{r/x)^f + (x, t)dx 

+ J {1 + s)^~^ j [{{v/x):^f + vl^ dxds < 

Step 3. In this step, we prove (13.1071) . Choose / = 1/2 in fl3.115p to get 



due to (13.41) . Then, (I3.107P follows from (I3.39P and (I3.116p . 


(3.114) 


(3.115) 


(3.116) 


(3.117) 


□ 


Lemma 3.9 Suppose that (13.ip and (13.21) hold. Let 6 G (0,7/2], and a and {3 be given 
respectively in (12.181) and (12.191) . If t E (0, (27 — 2 — 9)/%] and I E (0,1), then for t E [0,T], 


(1 + «)«-■ J (p«?) (x,t)dx + j\l + sf-' jf (vl 

(1 + tf-^ [ [{{v/x)^f + vl^] (x, f)dx < a,i^(0). 
Jo 


+ {vjxf) dxds < a^(0), (3.118) 

(3.119) 


Proof The proof consists of two steps. We prove (I3.118P and (13.1191) in Steps 1 and 2, 
respectively. 

Step 1. In this step, we prove that 


(1 PtY ^ J {'f’pVt) {x,t)dx + y (1 + s)^ ^ J 


(3.120) 


where ip is & smooth cut-off function defined on [0,1] satisfying (I3.110p . Then, (I3.118P follows 
from (^m> and (imH . 
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Differentiating fl3.6p with respect to t yields that 


p(-) vu-2p{^'\ \t-l 
\r J \r / X 


( 2 - + — 

r 




=iy 


r 


-l 

r / 2 2 \ 1 


/ \ 2 ' 

— z/ 

P' + 

-Ou 

p9 f^ + 2 -) 

- X 

L Vx r J\ 

X 

\rx r J 


+ Aui 6 


ft I Vx V 

pU— + 2- 
r 




(3.121) 


where p = p{r{x,t),t) = pr^^(r/x) Multiplying equation fl3.12ip by 'ijjvt, integrating the 
product with respect to the spatial variable, and noting that Vt{0,t) = v{0,t) = 0, one has 


d 

dt 


pij j v^dx + v f p^ (— + 2—'j {i/JVt)^ dx = ^ Jk, 
r J \rx r j 


(3.122) 


where 


~ J ~P^ {~^ v^dx + y v(()p'ipvtdx < C j p® dx, 


J2 = -l / P^ ( 2^ + ^ ) {:ipvt),dx, 


Js =1^ / p" 


^ f Vx 

+ 2— 0 [ -h 2 — 

\rx r 


i'lpvt)^ dx, 


J4 = - f p^(— +2-) Uvt-) dx, 
J \rx rj V r/x 


J 5 =4z/i / p^ 


ly^ j V 7^ fy' ^ ' 


dx. 


In view of (I3.95p . '^'{x) < 0 and (p^)x < 0, the second term on the left-hand side of 
(13.122p can be estimated as follows: 


P' (^ + 27) (^D),dx 



= / p^^jj—dx + 2 1 p^+vtVtxdx / p'^'i!)' ( — 2— ) vtdx 

J Tx Jr J \rx r J 

= /pV f — + - H{t) + [ pV f — + 2^) Vtdx, 

J \rx ry J \rx rJ 




6 I f ( o 


where 


and 


H{t) = j ti^y^dx + [(/)^ + Wx] vldx < J ^Wxv'fdx, 


Wx = 


P 


X 


2 \ 9 


j. 2 r 


- 1 


= ep\l + 0.f-^ilx + ^e{x,t). 


(3.123) 
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It follows from fl3.108p . fl3.9p . and the Holder and Cauchy-Schwarz inequalities that for any 

w > 0, 


H(t)< 


^^-v^dx + C [ iljWlvfdx 


.2 


1 

'4 

1 

'A 


^vldx + 2C 




tPWl 


+ C{ I ipwldx 


„2 



VtVtydy dx 


1 - 1/2 


<- 




»i-z 


X 2 


.2 


vfdx + Cioo / ip-^dx + Ciu 


< 


1 f 2 


X 



x^vl-^dx 


vfdx + Ci / p^{x^v^^ +v/)dx 


This, together with (I3.123p . implies 

j dx>c j ii (^vl^ + ^'^dx-Ci j / {x^vl^ + v^) dx. (3.124) 


Thus, by (I3.122p . one has 


d 

dt 




It therefore follows easily from the Cauchy-Schwarz inequality that 


d 

dt 



ipp' 


-e 


'xt 


+ -| ) dx 


<Ci 


+ vl + {v/xf + vl) dx. 


(3.125) 


So, fl3.120p follows from fl3.125p . fl3.63p and fl3.107p . 

Step 2. In a similar way to deriving (13.114p . one can use fl3.112p . fl3.120p . fl3.39p and 
fl3.63p to obtain 

(1 + t)^~^ j 'ip{x)p^'^~^{x)£if^{x,t)dx < C^S’{0). 

This, together with fl3.10ip and fl3.112p . gives fl3.119p . 

□ 


Lemma 3.10 Suppose that (13.ip and (13. 2 p hold. Let 9 G (0,7/2], and a and (3 he given 
respectively in (I2.18p and (I2.19p . If i & (0, (27 — 2 — 6 *)/ 8 ] and I G (0, 1), then for t G [0, T], 

(1+f)^"^ ||(xa:, v/x, r/x - l)(-,f)||5^i(jo;]) < C,,z(f(0). (3.126) 
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Proof. In a similar way as for fl3.117p . one can use fl3.39p . fl3.63p . fl3.106p and fl3.119p to 
obtain 


»l/2 


(1 + W [vl + {v/xf + {vx - 1)^ + {r/x - 1)^] {x, t)dx < C'.tf (0), (3.127) 


which implies, by using (I3.39p and (I3.63p again, that for any I G (0,1), 

(1 + f [vl + {v/xY + (tx - 1)^ + (r/x - 1)^] (x, t)dx < Q/^(0). (3.128) 


This, together with fl3.106p and fl3.119p . gives fl3.126p . 


□ 


3.3 Global existence of strong solutions 

In this subsection, we prove Theorem 12.21 by hrst verifying the a priori assumptions fl3.ip 
and fl3.2p . Choose i = (27 — 2 — d)/8 and / = 1/2 in fl3.126p to get 

IKu^;, v/x, rx-l, r/x- l)(-,t)||5^i([o^^/2]) ^ 

which implies, using || • H^oo < || • ||^i, that 

||(u^, v/x, r^-1, r/x-l)(-,f)||5^oo([o,i/2]) < C'^(O), tG[0,T]. 

It follows from choosing l = (27 — 2 — 6*)/8 in (I3.75p that 

||(u^, r/x, rx - 1, r/x - !)(•,i)|lioc([i/ 2 ,i]) < C^(0), t G [0,T]. 

This verihes the a priori assumption for small ^(0). 

The local existence of strong solutions for the problem fl2.12l) can be obtained easily 
by combining the approximation techniques in [T3] and the a priori estimates obtained in 
Sections 13.11 and 13.21 at least for small ^(0). Indeed, the a priori estimates obtained in 
Sections l3.1l and [T^ are also sufficient for the local existence theory, for small ^(0). Therefore, 
a standard continuation argument proves Theorem 12.21 with the estimates obtain in Sections 
IQ and 

Remark 3.11 The local existence theory in which is for the case of 9 = 1, does not apply 
to our case. 


4 Nonlinear Asymptotic Stability 


This section is devoted to proving Theorem 12.31 

First, it follows from the fact that ||5'|||oo < ||5'||L2||5'a:|lL2 for any function g, fl3.63p . fl3.64p . 
firm and fl3T5D that 


Y+tr 


„ r 37-2+2(n-e) a 1 
1 2 ( 7 + q :- 0 ) r 2 ’ 


7—1+ q : —6 
7+<a —6 


r(x,t) — xp < G^(0), X G /, 


(4.1) 
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(1 + t) 2 + x(l + t)^-40-4(7-l)-a} 
Since for small t, 


\u{r{x,t),t)\^ < ClS’{Q), X G I. (4.2) 


37-2 + 2(a-g) _ 1 ^ 7-1 + a-g 

2(7 + « — 6 *) 2 “ 7 + 0 — 0 


then we have 02.211) and 02.241) . 

Due to the fact that || • ||loo < || • and 03.126p . we have for any fixed I G (0,1), 

(1 + tf~^ {\vx{x,t)f + \x~^v{x,t)\‘^ + \ra:{x,t) - 1|^ + \x~^r{x,t) - 1|^) 

< C^S’{0), X G [0, /], 

which implies 02.23p . To prove 02.22p . we notice that for any 6 G [0, 2 — 7], 
x^p~^{x)\p{r{x,t),t) — p(x)p = x^p‘^~^{x)£f{x,t) 


= I iy^p^-^n\dy < 3 / y^f-^£l^dy + 2 / y^p'^-^aHydy 


2-2-bc^2, 


.3^2-6. 


<C / y‘^^{y^\ry - Ip + |r - y\‘^)dy 
Jo 

+ C( [ p^^-^Uldy] ( f yV-^'^-^^n^dy ) , 




\ 1/2 

V 


where the first inequality is due to 02.31) : and if 7 > (5 — 26)/3 

J y^p^~^’'~^^£l^dy < (^J - Ip + \r- y\^)dy 

X f/ 2/V(2/^|r,-lp + |r-|/p)dJ 


5-27-0-26 

7-0 


37-5+26 

7-0 


Then, it follows from 03.39p and 03.1O8P that for any be [0, 2 — 7], 

o L /3 max{0, 37—5+26} 

x^p“^(x)(l + f) 2 2(7-s) |p(r(a:, f), t) — p(x)P < C<f (0), xG[0,1]. 

Due to 04.31) and 6 < 1, we have 

{l + tf~^p~’'{x)\p{r{x,t),t) -p(x)p < C^S’{0), X G [ 0 ,/]. 


So, 02.221) is a consequence of 04.4p and 04.51) . 

It remains to proving 02.251) . We use 03.92p . 03.84p . 03.88l) - 03.90p . 03.93P and 

1/2 



(4.3) 


(4.4) 


(4.5) 
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to get 


‘^\x)\r^{x,t) -1\^ + C / y'^p'^ ^vl{y,t)dy 


+ Ct 


(1 + + (1 + + «- 4 (T-l)-o)j gf^y 


(4.6) 


Here (13.631) . fl3.64p and fl3.39p have been used. In view of fl3.64p . we see that if 7 — 1 > 0, 

J y^p'^~^vl{y, t)dy < C,{1 + ( 0 ); 

and if 7 — 1 < 6 *, 

[ y‘^p'~^vl{y,t)dy <C { I y^p%l{y,t)dy 


a-29+2(7-l) 28-2(7-1) 


< 0(1 + 

That means 

/ <0(1 + 2»-2h-i)).f(0), 


(4.7) 


In a similar way to deriving (14.4p . we have 

x^p^'^-^\x)\r^{x,t) - 1|^ < C,(l = (4.8) 

So, (12.251) is a conclusion of (I4.6l) - (l4.8p and (14.3p . 
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